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Abstract 

The production of two high-p^ jets in the interactions of quasi-real photons 
in e + e~ collisions at -/s^ from 189 GeV to 209 GeV is studied with data 
corresponding to an integrated e + e~ luminosity of 550 pb -1 . The jets re- 
constructed by the fcj_-cluster algorithm are defined within the pseudo-rapidity 
range — 1 < i] < 1 and with jet transverse momentum, p?, above 3 GeV/c. The 
differential di-jet cross-section is measured as a function of the mean transverse 
momentum p? of the jets and is compared to perturbative QCD calculations. 
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1 Introduction 

This paper presents the results of a study of di-jet events produced in two-photon 
collisions in the anti-tagged mode, i.e. when both scattered electrons 1 escape detection. 
Large px processes involving quasi-real photons are sensitive to both, quark and gluon, 
components of the resolved photons. Thus, the analysis of high-p T jet production com- 
plements the studies of the deep-inelastic structure function of a quasi-real photon which 
probe the quark distribution. Considered together, they allow the parton density function 
of the photon to be determined. The perturbative QCD scale of the hard interactions is 
provided by the jet transverse energy to be used in the calculations. Available leading and 
next-to-leading order calculations [1] can be tested with the large samples of LEP-2 data. 
The present analysis adds new results to those obtained by other LEP experiments [2,3]. 

Gamma-gamma collisions exhibit a behavior typical of hadron-hadron interactions, i.e. 
the centre-of-mass frame of the hard scattered partons is moving in the 77 centre-of-mass 
frame. The /cj_-cluster algorithm [4-6], invariant under boosts along the beam axis, is 
used. The variables adopted for the jet analysis are the transverse momentum pt, the 
azimuthal angle <fi and the pseudo-rapidity rj = — m(tan((9/2)) of the particle(jet) 2 . 

Hadron production in the collisions of quasi-real photons is described by the set of 
processes illustrated in Fig. 1. The interaction of bare photons (the direct term) is 
described by the Born-box diagram within the quark-parton model (QPM). If one or both 
photons are resolved into a partonic structure, the process is termed single- or double- 
resolved, respectively. A part of the double-resolved interactions with both photons 
resolved into a bound quark-pair system is described by the vector dominance model 
(VDM). The relative contribution of the different components depends on the kinematic 
regime. The present study involves such a 3-component model in order to describe the 
data and to correct on the experimental effects and then to compare the data with leading 
order (LO) and next-to-leading order (NLO) QCD calculations. 

An important issue is the expected difference in topology. Almost all hadrons pro- 
duced in QPM-like events should belong to two hard-jets, while two jets in events with 
resolved photons are accompanied by remnant jet(s). The variables sensitive to such a 
difference are [7] 

X)(£jet+Pz,jet) ZX-Ejet-PzJet) 

^,-\- j c ts ctnd oc — j°^ s 

7 Yl (£part+£>z ; part) ' ' ' 7 J2 (-^part "pz,part) ' 

part part 

where 'part' corresponds to all detected particles and E- set and p z ,\ct are the two hard-jets 
energy and the component of jet momentum along the z-axis, respectively. These variables 
represent the respective fractions of the 7 momenta relevant to the hard interaction. The 
photons in the QPM-like events participate in the interaction entirely and both x+ and 
x~ should be equal to unity, while the presence of a remnant jet (single-resolved photon) 
shifts x+ or x~ toward some lower values. In the double-resolved case, both variables 
are far from unity. 

Experimental details are discussed in Section 2. Section 3 presents the analysis of the 
data which are compared to the leading and next-to-leading order QCD calculations [1] . 
Results and conclusions are given in Section 4. 

1 Throughout this paper, electron stands for electron and positron. 

2 The origin of the DELPHI reference system is at the centre of the detector. The z-axis is parallel to the e~ beam, the 
x-axis points horizontally to the centre of the LEP ring and the y-axis points vertically upwards. The coordinates R, <f>,z 
form a cylindrical coordinate system and 6 is the polar angle with respect to the z-axis. 
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2 Detector and Data 

The DELPHI apparatus and performance are described in detail in ref. [8,9]. 

Charged particles are detected in the tracking system containing a 3-layer silicon Ver- 
tex Detector, the Inner Detector, the Time Projection Chamber and the Outer Detector. 
The Forward Chambers extend the acceptance into the forward region. The tracking 
detectors are located inside the superconducting solenoid providing a uniform magnetic 
field of 1.23 T parallel to the axis of the colliding e + /e~ beams. The combined momen- 
tum resolution provided by the tracking system is of the order of a few per-mille in the 
momentum range of the present study. 

Photons are detected in the electromagnetic calorimeters - High density Projection 
Chamber in the barrel region and the Forward Electro-Magnetic Calorimeter in the end- 
caps with polar angle coverages of (40° — 140°) and (9° — 35°, 145° — 171°), respectively. 

Jets are reconstructed from two kinds of detected objects - charged particles and 
photons. 

The events of interest are quite energetic and the multiplicity of final state particles is 
large. The events are triggered by several components of the DELPHI trigger system [10] 
and the resulting efficiency is estimated to be equal to unity. 

The data used for this analysis have been collected by the DELPHI detector at LEP 
at an e + e~ centre-of-mass energy (y/s^) from 189 GeV to 209 GeV and correspond to 
an integrated e + e~ luminosity of 550 pb _1 . 

2.1 Data selection 

The following criteria are applied to select the data sample. 

• Multihadron event selection. 

Each event should contain at least 5 charged particles. A charged particle is counted 
if its momentum is greater than 0.2 GeV/c, the polar angle is within the interval 
from 20° to 160°, the relative momentum error is less than 100% and the impact 
parameters are smaller than 4 cm in Rip and 10 cm in z. 

• Background suppression. 

The invariant mass of the system calculated from the charged particles (assuming 
them to be pions) and the photons (detected in the electromagnetic calorimeters) 
should be below 35 GeV/c 2 . The energy threshold for an electromagnetic calorimeter 
cluster is set to 0.5 GeV. The total transverse momentum | XX~Pr)| should be below 
30 GeV/c. These criteria suppress events coming from e + e~ annihilation and they 
set upper limits (see below) for jet transverse momentum. 

• Anti-tagging condition. 

There should be no clusters in the luminosity monitor STIC (which covers the region 
from 29 mrad to 185 mrad in the polar angle) with energy greater than 25 GeV (anti- 
tagging condition). This ensures the photons to be quasi-real. 

• Di-jet event selection. 

Jets are reconstructed by the /c^-cluster algorithm, implemented in the KTCLUS 
program [6]. It operates with objects (at generator level they may be either particles 
or partons) and tries to unite them into a jet in a recursive way. The following 
variables are calculated for each object i and a pair of objects i and j: 

di = E\ { and d {j = min(d { , d^lfa - rjj) 2 + (fa - (p^/R 2 , 

where the transverse energy Et,% is calculated with pion mass assumption for a 
charged particle and zero mass for a photon, R is usually set to unity. If the smallest 
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value (d min ) of any d in this sample is the dij value of a pair then these two 

objects are joined into a new object and its E T , <p and rj are recalculated. Otherwise, 
if the smallest value in the sample is the di value of an object i, the corresponding 
object is not mergeable and it is moved from the object list to the jet list. The 
procedure continues until there are no more objects. The pseudo-rapidity rj of exactly 
two jets must be in the range — 1 < rj < 1 and their transverse momenta, px, above 
3 GeV/c. Besides these two hard process jets, the event may contain one or more 
jets outside the pseudo-rapidity and px domain mentioned above. They are termed 
'remnant' jets. The choice of the kinematic limits is conditioned by the desire to 
suppress soft 77 interactions and to keep hard process jets within the DELPHI 
acceptance. The cut on jet pr should not be too large since, besides decreasing the 
measurement accuracy, it reduces the contribution coming from events with resolved 
photons. The extension of the rj domain results in an increase of threshold effects, 
since some of the produced jet particles are outside the detector acceptance and the 
efficiency of the jet reconstruction is decreased. 

An additional cut on the mean jets transverse momentum pr > 4 GeV/c allows a 
Pt asymmetry of the reconstructed jets. If such a cut is applied symmetrically to 
each of the two jets, the QCD calculations of some variables are unstable [1]. 

The above criteria select a data sample of 5147 events. The data are corrected for 
detector inefficiency and acceptance effects with a leading-order (LO) Monte Carlo simu- 
lation described below. LO simulation is also used for the estimation of the soft processes 
accompanying the hard-initiated jets and for taking into account the influence of hadro- 
nisation. 

2.2 Monte Carlo simulation and background estimation 

The Monte Carlo generator PYTHIA (version 6.205) [11] is used for the simulation of 
hadron production in 77 interactions. The program contains an interface to an external 
library of parton density functions (PDF) for the photons [12]. The default SaSID [13] 
PDF is taken. The soft underlying events are modeled through multiple parton interac- 
tions (MIA) of several parton pairs within the same event. An event with both photons 
resolved may contain, besides hard initiated jets and the corresponding remnant jets, 
such an additional MIA contribution. 

The simulated sample is approximately 4 times larger than the data sample. 

The main background process (e + e~ — > hadrons) is simulated by the KK2f event 
generator (version 4.14) [14]. The expected number of these background events in the 
selected data sample is estimated to (500 ± 5) events. The contamination of r pairs 
produced in two-photon interactions is evaluated as (43 ± 3) events, using the BDKRC 
program [15]. The background of r pairs produced in e + e" annihilation is negligible. The 
events with W, Z bosons contribute (38 ± 4) events. 

Note that the mentioned background sources (they will be referred to as 'non-77' 
background) are only a part of the total background. Additional sources will be discussed 
later in Section 3.2. 

3 Analysis 

The entire x~)-space can be split into four quadrants (Fig. 2) by x 1 equal to 0.85. 
The chosen value of x 7 selects approximately equal statistics in the different domains. 
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One can define three kinematic regions: the first - x^ and x~ greater than 0.85 (the 
region is termed as 'Dir' domain), the second - when both values are below 0.85 (double- 
resolved domain, 'DR') and the third region - when one variable is below 0.85 while the 
other is above it (single- resolved, 'SR'). 

Figure 3 shows the total energy outside the reconstructed jets (E out ) in comparison 
with the model expectation. The contributions of different parts of the model are also 
presented in the figure. 

The 'Dir' domain (Fig. 3b) is mostly populated by QPM-like events (94%), the con- 
tribution of the double- resolved subprocesses to the 'DR' domain (Fig. 3d) is around 
54%, while the 'SR' domain (Fig. 3c) contains all three types in nearly equal parts. It is 
seen that the model does not describe the data in the parts of the x~)-space where 
the contribution of the resolved processes is essential. It is also clear that a simple one- 
parameter renormalization of the model as a whole as done in previous publications [2,3] 
is not adequate since its components have to be tuned separately with different factors. 
This renormalization is a new approach to the analysis of jets produced in 77 collisions. 

A 3-parameters fit of the data is performed simultaneously on the following distribu- 
tions: E out (Fig. 3), the transverse momentum balance ( IX) Pi^l) ( n °t shown) and 
the total invariant mass W of the detected particles (not shown). The distributions are 
fitted in each x~)-domain. The obtained parameters (scale factors a for each model 
component) are a qpm = (0.86 ± 0.02), a s - res = (1.49 ± 0.09) and ad-res = (1-93 ± 0.05). 

Figure 4 presents the background-subtracted data compared to the 77 simulation 
(with scaling factors applied) for the variables which were not involved in the fit. The 
data are in good agreement with the simulation. The Monte Carlo set, scaled according to 
the fit parameters, consists of 33% QPM-like events, 23% single-resolved and 44% double- 
resolved events. Double-resolved events are subdivided into contributions 'no-MIA' (75%) 
and 'with-MIA' (25%). 

3.1 Multiple parton interactions 

As already mentioned the modeling of the double-resolved events includes the so-called 
multiple parton interactions (MIA) [11]. Since the resolved photons are composite objects 
consisting of many partons, it is assumed that different pairwise interactions may take 
place during one 77 collision. It has been shown that the inclusion of MIA improves 
the description of the data [16]. The amount of MIA in PYTHIA is determined by a 
cutoff parameter on the parton transverse momentum. The default value of 1.6 GeV/c 
is used. MIA influence the model predictions in two ways. The first is the case when the 
hard initiated jets satisfy the selection criteria and MIA systematically increase the jet 
transverse energy. The second way is when the hard process alone does not provide two 
jets passing the selection criteria but combined with MIA initiated particles the event 
is selected. The MIA contribution in this case is to the background and it has to be 
estimated and subtracted from the data according to the model expectation. For this 
purpose PYTHIA is run in the mode with MIA switched off. It is found that 27% of 
the double-resolved 77 interactions with MIA contributions would not pass the di-jet 
selection criteria if MIA were absent. Note that the double-resolved part of the model is 
fitted in a coherent way, i.e. without splitting it into subsamples with-MIA and without 
it. 

Summarizing the MIA influence, the model predicts the background coming from MIA 
as 7% (3%) of the double-resolved part (the total model prediction). On the other hand 
18% of the double-resolved events are affected by MIA via the increase of the jet transverse 



■5 



momentum. The absolute value of the background coming from the MIA contribution is 
estimated as (131 ± 7) events. The importance of the 3-parameter fit performed has to 
be emphasized since the evaluation of the background is done by using the scale factors 
obtained. 

3.2 Hadronisation corrections, acceptance and background sub- 
traction 

The theoretical predictions have to be transformed from the parton level to the level at 
which particles are produced in the 77 interaction vertex - the particle level ('hadronisa- 
tion corrections'). The Monte Carlo partons are considered in leading-order and are not 
identical to the NLO partons in the theoretical calculations. The corrections have been 
calculated using the PYTHIA program. The distribution of the mean transverse momen- 
tum of the jets pr (difference of jet pseudo-rapidities | A77I) obtained with the /cj_-cluster 
algorithm at the parton level, is divided by the distribution obtained after hadronisation. 
The resulting 'hadronisation corrections' are shown in Fig. 5. They reach 60% value at 
smallest pt and decrease to 30% for higher transverse momentum. 

The next step is to transform the data from the level of the detected particles to 
the hadrons produced in 77 interactions in order to take the detector acceptance into 
account. The detector acceptance enters in two ways - either some events do not satisfy 
the selection and thus they are lost or the observed distributions in the selected events are 
distorted. Bin-by-bin corrections are used. Figure 6 illustrates the detector efficiency as 
a function of pr (the dip in the middle of the distribution is caused by limited statistics 
of the simulation) and 

Note that these corrections apply to two-jet events. The first correction transfers the 
two-jet calculations carried out at the parton level to model expectations at the level 
of the produced particles. The acceptance has been calculated both at the produced 
and detected particle levels for the two-jet events. That is why the following and last 
source of background ('non2-to-2 jets') should be taken into account. Some events are 
reconstructed with two jets because undetected particles change the event topology but 
they do not have two jets at the production level. This kind of background is estimated 
as (893 ±13) events and it becomes the largest source of background for the di-jet studies. 

4 Results and conclusions 

Figures 7(a-b) show the pr and \Arj\ data distributions compared to the total simu- 
lation. All sources of the background mentioned above are given there as well. 

The data are bin-to-bin corrected for the acceptance and, applying the integrated lu- 
minosity of LEP corresponding to the data, the total di-jet cross-section is measured to 
be 17.8 ± 0.6 pb for jets within the pseudo-rapidity range of — 1 < rj < 1 and for a jet 
transverse momentum p T above 3 GeV/c. The statistical and systematic uncertainties 
(approximately of equal size) are added in quadrature. The systematic uncertainty in- 
cludes the errors of the fit mentioned above, the uncertainty coming from MIA handling 
and the background estimations. The model expectations are 20.5 ±0.1 pb and 18.1 ±0.6 
pb for the calculations carried out in the leading and next-to-leading order, respectively. 
The differential cross-sections are obtained with bin-by-bin corrections for the detector 
efficiency (Fig. 6) for the data and with 'hadronisation corrections' (Fig. 5) applied to the 
theoretical calculations. The obtained cross-sections are shown in Fig. 8 and compared to 
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the predictions mentioned above (the NLO calculations done with Monte Carlo method 
have sizable uncertainties which are shown on the plots). Numerical results are presented 
in Table 1 as a function of pr and in Table 2 as a function of |A^|. The errors include 
statistical and systematic uncertainties added in quadrature. 

Since we scaled the contribution of QPM-like events ('direct' term) and this part is 
well calculable, it seems interesting to compare the scaled distributions to the results 
of calculations carried out in the leading and next-to-leading orders. The results are 
presented in Fig. 9. The scaled 'direct' term is only around 10% below NLO calculations; 
the main difference comes from the small pr domain, while there is good agreement with 
NLO calculations for pff greater than 4.5 GeV/c. 

In conclusion, the production of two high-px jets in the interactions of quasi-real 
photons is studied with the DELPHI data taken at LEP-2 with an integrated e + e~ 
luminosity of 550 pb _1 . The jets reconstructed by the fcj_-cluster algorithm are defined 
within the pseudo-rapidity range of —1 < rj < 1 and for a jet transverse momentum p? 
above 3 GeV/c. The total and differential di-jet cross-sections are measured for a mean jet 
momentum pr from 4 GeV/c to 14 GeV/c. The total cross-section agrees with the next- 
to-leading order perturbative QCD calculation within the experimental uncertainties and 
is 18% below the calculation carried out in the leading order. The measured differential 
di-jet cross-section is also found to be in good agreement with NLO QCD predictions [1] 
and complements results obtained by other LEP experiments [2,3]. 
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VT ,mean 


&rneasured,total (pb/GeV) 


a LO (pb/GeV) 


onlo (pb/GeV) 


4.0-4.5 


9.59 ±0.53 


11.59 ±0.13 


9.75 ±0.85 


4.5-5.0 


6.71 ±0.41 


7.61 ±0.10 


6.71 ±0.53 


5.0-5.5 


4.56 ±0.33 


5.24 ±0.07 


4.49 ±0.27 


5.5-6.0 


3.06 ±0.26 


3.86 ±0.06 


3.49 ±0.19 


6.0-7.0 


2.19 ±0.16 


2.46 ± 0.04 


2.28 ±0.12 


7.0-8.0 


1.17 ±0.12 


1.40 ±0.03 


1.28 ±0.08 


8.0-9.0 


0.88 ±0.12 


0.88 ±0.02 


0.79 ±0.06 


9.0-10.0 


0.71 ±0.11 


0.61 ±0.02 


0.58 ±0.05 


10.0-12.0 


0.31 ±0.05 


0.33 ±0.01 


0.31 ±0.03 


12.0-14.0 


0.11 ±0.03 


0.18 ±0.01 


0.16 ±0.03 



Table 1: Measured cross-section of di-jet production in quasi-real 77 interactions as a 
function of pr together with leading order and next-to- leading order calculations [1]. 



|Ai7| 


& measured, total (pb) 


o"lo (pb) 


onlo (pb) 


0.-0.2 


3.50 ±0.21 


4.23 ±0.05 


3.44 ±0.25 


0.2-0.4 


2.96 ±0.18 


3.66 ±0.05 


2.76 ±0.20 


0.4-0.6 


2.54 ±0.18 


3.31 ±0.04 


3.32 ±0.25 


0.6-0.8 


2.67 ±0.19 


2.81 ±0.04 


2.36 ±0.19 


0.8-1.0 


2.09 ±0.17 


2.35 ±0.04 


2.20 ±0.24 


1.0-1.2 


1.78 ±0.16 


1.87 ±0.03 


1.84 ±0.20 


1.2-1.4 


1.39 ±0.14 


1.38 ±0.03 


1.28 ±0.13 


1.4-1.6 


0.67 ±0.11 


0.95 ±0.02 


0.80 ±0.09 


1.6-1.8 


0.29 ±0.12 


0.54 ±0.02 


0.63 ±0.09 


1.8-2.0 


0.29 ±0.11 


0.16 ±0.02 


0.09 ±0.04 



Table 2: Measured cross-section of di-jet production in quasi-real 77 interactions as a 
function of |A^| together with leading order and next-to-leading order calculations [1]. 
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Figure 1: Main diagrams corresponding to the hadron production in 77 interactions via 
vector meson interactions (VDM-like, a), point-like interactions (QPM-like, b) and with 
one (c) or both (d) photons resolved into partons. 
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Figure 2: (x+, )-space divided by x 7 into four domains ('Dir', 'SR' and 'DR' notations 
are introduced in section 3. 
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Figure 3: Comparison of the total energy outside the reconstructed jets E out with the 
simulation for the selected data (a), and data in the phase-space domains 'Dir' (b), £ SR' 
(c) and 'DR' (d), described in section 3 of the text. The thick solid histograms correspond 
to the sum of the di-jet 77 interactions as predicted by PYTHIA and the background 
processes. The contributions of different 77 subprocesses are shown by dotted (QPM 
term), dashed-dotted (single- resolved) and thin solid (double-resolved) histograms. 
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Figure 4: Comparison of the background subtracted data with the simulation after the fit 
for the charged multiplicity (a), the invariant mass calculated with the charged particles 
only (b), (x+ + x~)/2 (c) and the mean number of particles in the reconstructed jets (d). 
The simulation (thick solid histograms) sums all three model components scaled with 
the factors given in the text. The contributions of 77 subprocesses are shown by dotted 
(QPMterm), dashed-dotted (single-resolved) and thin solid (double-resolved) histograms. 
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Figure 5: Hadronisation corrections as a function of the mean jet transverse momentum 
Pr (a) and the absolute difference of jet rapidities \ Ar)\ (b). 
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Figure 6: Detector efficiency as a function of pr (a) and |A^| (b). 



14 



DELPHI 



DELPHI 



^2500 




s 

&1000 



800 - 



600 



400 



200 - 



b) 



" ' ■ ■ - - . X- l-l-l J 



10 11 12 13 14 
Pr,,nean(GeV/c) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 7: Distributions of the mean jet transverse momentum (a) and the absolute 
difference of jet rapidities \ Ar]\ (b). The data are presented by bars, the sum of the scaled 
PYTHIA prediction and of the background by solid histograms. The background com- 
ponents are shown by dash-dotted lines (MIA contribution), by dashed lines ('non2-to-2 
jets' background) and by shaded histograms ('non-77' background). 
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Figure 8: Cross-section of di-jet production in quasi-real 77 interactions as a function 
of (a) and |A^| (b). The jets are reconstructed by the /cj_-cluster algorithm within 
the pseudo-rapidity range of —1 < rj < 1 and the jet transverse momentum above 3 
GeV/c. Dashed histograms show the leading order calculations and next-to-leading order 
calculations [1] are presented by open diamonds. 
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Figure 9: Cross-section of the 'direct' term as a function of pr (a) and |A^| (b). Dashed 
histograms correspond to the PYTHIA predictions scaled according to the fit. The circles 
(squares) present the results of the NLO (LO) calculations [1]. 



